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ABSTRACT 


A method is described for obtaining optimal attitude estimation algorithms for 
spacecraft lacking attitude rate measurement devices (rate gyros), and then demonstrated 
using actual flight data from the Solar, Anomalous, and Magnetospheric Particle Explorer 
(SAMPEX) spacecraft. SAMPEX does not have on-board rate sensing, and relies on sun 
sensors and a three-axis magnetometer for attitude determination. Problems arise since 
typical attitude estimation is accomplished by filtering measurements of both attitude 
and attitude rates. Rates are nearly always sampled much more densely than are 
attitudes. Thus, the absence/loss of rate data normally reduces both the total amount 
of data available and the sampling density (in time) by a substantial fraction. As a 
result, the sensitivity of the estimates to model uncertainty and to measurement noise 
increases. In order to maintain accuracy in the attitude estimates, there is increased need 
for accurate models of the rotational dynamics. The proposed approach is based on the 
Minimum Model Error (MME) optimal estimation strategy, which has been successfully 
applied to estimation of poorly-modeled dynamic systems which are relatively sparsely 
and/or noisily measured. The MME estimates may be used to construct accurate models 
of the system dynamics (i.e. perform system model identification). Thus, an MME- 

based approach directly addresses the problems created by absence of attitude rate 
measurements. 


t 


% 

9 


NASA Graduate Research Fellow; Student Member AIAA 
Department of Mechanical and Aerospace Engineering, State University of New York 
NRC Postdoctoral Research Fellow; Member AIAA 
Goddard Space Flight Center, Code 712, Greenbelt, MD 20771 
Postdoctoral Research Assistant; Member AIAA 
Department of Mechanical and Aerospace Engineering, State University of New York 
Research Assistant; Student Member AIAA 
Department of Mechanical and Aerospace Engineering, State University of New York 
Associate Professor; Senior Member AIAA 
Department of Mechanical and Aerospace Engineering, State University of New York 


at Buffalo, Buffalo, NY 14260 


at Buffalo, Buffalo, NY 14260 
at Buffalo, Buffalo, NY 14260 
at Buffalo, Buffalo, NY 14260 


PR6CWWNG PAGE BLANK NOT FH.MRO 


497 


INTRODUCTION 

The attitude of a spacecraft can be estimated by either single-frame deterministic 
methods (such as TRIAD and QUEST [1-2]) or algorithms which combine analytical 
models with attitude measurements and, for most spacecraft, attitude rates (such as the 
Kalman filter [3]). Generally, the use of rate gyros significantly improves the attitude 
estimation, because the densely-measured rates may virtually eliminate the need for 
dynamic models. However, the intentional omission of rate gyros in the design of 
satellites is increasingly likely as resources become more scarce (for example, the Solar, 
Anomalous, and Magnetospheric Particle Explorer (SAMPEX) satellite does not have rate 
gyros on board). In addition, existing satellites with rate gyros on board may experience 
gyro degradation or failure (such as the failure of four of the six rate gyros on the Earth 
Radiation Budget Satellite (ERBS) [4]). When rate gyro data is either not available or not 
dependable, the attitude estimation accuracy becomes much more heavily dependent on 
accurate dynamic models since the attitude measurements are typically much less dense 
and less accurate than the rate data. In these cases, dynamic models may be required to 
provide estimates between and/or in addition to the attitude measurements. 

Unfortunately, accurate models of spacecraft rotational dynamics are often unavail- 
able. In cases where the spacecraft was launched with rate gyros, the attitude estimation 
algorithm likely did not require an accurate dynamics model since dense rate measure- 
ments were available. Even for spacecraft which do not have rate gyros, determining an 
accurate rotational dynamics model may be difficult. If an accurate model is necessary in 
the attitude estimation algorithm, estimation accuracy is compromised. This is especially 
true for spacecraft launched with rate gyros which subsequently fail. 

To circumvent the problem of rate gyro omission or failure, analytical models of gyro 
biases can be used. An example of a commonly used gyro bias model is the model based 
on a Markov (exponential decay) process. This simple model has been successfully used 
in a Real-Time Sequential Filter (RTSF) algorithm in order to propagate dynamic state 
estimates and error covariances for the SAMPEX satellite (see [5]). A clear advantage 
to using dynamic models for gyro biases was shown for the case of Sun-magnetic near 
co-alignment. For this case, the single-frame algorithms, TRIAD and QUEST, showed 
anomalous behaviors with extreme deviations in attitude estimations. However, since the 
RTSF propagates an analytical model of the gyro bias, the attitude estimates are improved 
even when data from only one sensor is available (i.e., only magnetometer measurements). 

In theory, perfect, solvable models of the spacecraft rotational dynamics could be used 
to obtain perfect attitude estimates. When accurate rate gyros are present, they can often 
take the place of the dynamic models. When rate gyros are either absent or excessively 
noisy attitude estimation accuracy becomes critically dependent on the accuracy of the 
rotational dynamic models. The ERBS studies [6-8] showed, for an existing satellite, 
that modeling of the attitude dynamics leads to accurate attitude estimation algorithms. 
However, the authors concluded that in order to be operationally useful, “automatic” 
methods for determining these dynamic models must be available. 
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In this paper, a technique is described which directly addresses the problem of attitude 
estimation when rate data is not available (or severely degraded), regaidless of the cause 
e method described herein addresses this problem directly in two distinct but related 
approaches. First, the MME [9-12] method may be used simply to obtain accuram s^ 
estimates for dynaniic systems which are both poorly modeled and sparsely measured 

. A “° mp !' shed ,hrough ex P llclt accounting for errors in the dynamic model Thus 
nude estimation using existing satellite dynamic models (which may not be particularly 
accurate) is poss.ble. However, the MME estimates may also be used to consTrucUccurare 
models of the system dynamics (i.e„ perform system model identification). Urns the 
econd, and main thrust of the approach is the use of the MME to create more accurate 
ynamic models for use in ANY estimation algorithm (batch, sequential, or MME) 

An optimal attitude estimation algorithm is described which is capable of robust and 
“ “a', ° f SPa “f kin8 b °‘ h ™ rate Measurements 1 and 

Sd n he cumnt appr ° ach is based the Mi" 1 ®™ 

esrimr, c r M P esl,maUon “^gy. which has been successfully applied to 

estimation of numerous poorly-modeled dynamic systems which are relatively sparsely 

and/or noisily measured. The MME-based approach described in this paper has the 
capa 1 ity to automatically determine accurate rotational dynamic models resulting in 

as ZwTiT 191 e t, b “ ‘g 0 h ' 8h aCC “ raCy ° f esdmation “Stag accurate dynamic models, 
as shown m [9], while eliminating the practical limitations currently imposed bv the 

requirement that the models be determined manually for each orbit Y 

The organization of this paper is as follows. First, a brief description of the SAMPFX 
“ “^nations of morion is shot*,. Then, ^rief tm^ 

es imation algorithm for nonlinear systems is shown. An MME estimator 
which incorporates the SAMPEX model, is next applied to estimate The d “ s 
udes, angular rates, and angular momentum) of SAMPEX using actual telemetry 

~mT- astIy> candidate functionai foms for ^ 

y E estimator are investigated. Results are compared with actual telemetrydata. 


SAMPEX MISSION DESCRIPTION 

The Goddard Space Flight Center (GSFC) Small Explorer (SMEX) program was 
developed to provide relatively inexpensive, frequent space science missions The Solar 

SMEX T’ Magne i os P henc Particle Explorer (SAMPEX) satellite is the first of the 
SMEX class missions. The SAMPEX [5] general mission is to study energetic particles 

The spacecraft is three-axis stabilized in a 550 by 675 km elliptical orbit with an 82° 
inclination. The nominal mode is a one rotation per orbit about die Sun vector The bSy 
z-axts is defined by the instrument boresights and is required to be within 15° of zenith 
near the magnetic pole. The body y-axis nominally is aligned with the Sun vector 
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TTk attitude determination haidware consists of five coarse Sun sensors (CSS) 
(primarily ^or Sun-acquisition), one fine Sun sensor (FSS) (for roll and yaw), and a 
three axil magnetomerer (TAM) (for pitch). The attitude control hardware consists of 
a magnetic Mquer assembly (MTA) (for roll and yaw), and a reaction wheel assem y 
(RWA) (for pitch). The nominal attitude determination accuracy is 


SRTADS 

The SAMPEX Real-Time-Attitude-Determination-System (SRTADS) is a graphical 
uscrdnrerface program which compares and displays attitude solutions along with releme- 
« m^uremen™ in a real-time mode during real-time contacts. One of the functions 
J the system is to serve as a resting platform for filtering methods used in attitude 

determination. 

The current version of the SRTADS program implements three different attitud 

. C t . . n thp triad algorithm; 2) the QUEST algorithm; and 3) 

determination methods: 1) the ikiau aiguuuim, ; ^ , , 

the Real-Time Sequential Filter (RTSF). Both TRIAD and QUEST are single frame 
deterministic methods which primarily rely on a pair of measured vectors for attitud 
determinaUon. The RTSF is an extended Kalman filter which combines both measured 
data and a system model to obtain an attitude solution. 

THE MME APPROACH 

The Minimum Model Error (MME) estimation algorithm was developed for -optimal 
state estimation of poorly modeled dynamic systems ([9]). Motivated by problems 
satellite orbit/attitude determination (see [13]), in which significant unmo*M^namres 
may be present, the MME was formulated to rigorously account for both sigmfica 
modeling error and significant measurement noise. 

The MME state and model error estimations have been shown to be extreme y 
accurate in previous work [9-10], [12] and the algorithm shown to be f obust wtth 
respect to modeling errors, measurement errors, and measurement sparsity [12]. 
mediate trajectories are accurately estimated, *(() » JL«. and, most importan for 
the 6 realization/identification preblem, d(t) approaches the correct model error trajectory. 
Another key feature of the MME (explained in [101) is that the stare estimates are 
of jump discontinuities evident in Kalman filters, for example. 

The MME solution yields the optimal state estimates x(t) and the optimal model error 
estimates d(t ) Results presented in [12], [13] showed that for poorly modeled systems, 
the MME state estimates are of considerably higher accuracy than those obtained using 
standard^ approaches 1 based on Kalman filtering. In addition, the MME has been used 
as the basis for highly accurate and robust system identification algorithms, both linear 
[11], [14-16], and nonlinear [12], [17-18], based on the combination of state and model 

error estimates. 
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MODEL EQUATIONS 


The following is a brief summary of the kinematic and dynamic equations of 
motion for a three-axis stabilized spacecraft. The rotational orientation of the spacecraft 
(kinematic equations) may be represented by the quaternion attitude parameterization as 


£ = 

where 
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The elements of Q are the components of the instantaneous spacecraft angular velocity 
defined relative to the body frame. 

The dynamic equations of motion (Euler’s Equations of Motion) for a non-rigid 
spacecraft (SAMPEX is not modeled as a rigid body because it contains a reaction-wheel 
assembly), may be defined as 


dh 

dt 


N — u> x L 


— intody T" \Itw X W rw ]j 


(8a) 

(8b) 


where 

N = total external torque 
w = instantaneous angular velocity 
L = total angular momentum 
—body = angular momentum of the body 
drw inertia of the reaction wheel 
<jj TW = angular velocity of the reaction wheel 

Here, again, all vectors are resolved in a body-fixed coordinate system. The angular 
momentum of the reaction wheel only acts along the y-body axis. The nonlinear state- 
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space representation of the dynamic equations of motion is given by 
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The body angular rates can be determined using the angular momentum, reaction wheel 
momentum, and known inertia: 



Lx 

= I~ l { L 2 - Itw X w r 

L 3 


( 10 ) 


where I represents the inertia tensor ol the satellite. 

SRTADS provides time histories for all external torques (aerodynamic, gravitational, 
solar, magnetic, etc) and the reaction wheel input. The nominal numerical values for the 
SAMPEX inertia tensor and reaction wheel inertia are given by [5] 


I = 


'15.516 

0 

0 


0 

21.621 

-0.194 


0 

-0.194 

15.234 


kg - m 2 


I TW = 0.0041488 kg - m 2 


RESULTS 

A nominal satellite pass is used to compare the MME estimator to the deterministic 
method (e.g., TRIAD). Nominally, both FSS and TAM data is available throughout the 
orbit. Anomalous behavior occurs during either sun occultation and/or co-alignment of 
the measurement vectors. However, this test case involves a non-event pass (i.e., no 
anomalous behavior). SRTADS utilizes telemetry and ephemeris data from an orbital 
pass and determines the Euler attitudes using TRIAD, QUEST, and the RTSF in 5 
second intervals. 

The MME estimator uses a priori values from the single-frame solutions (i.e., the 
TRIAD solution for the spacecraft’s attitude). The MME estimator is then used to obtain 
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estimates for both the three-axis attitudes and the dynamic rates. Since this test case 
involves a non-event pass, the TRIAD solution for the spacecraft’s attitude is assumed 

qa unrv t Ct r (n0minaI) soIutlon - Figures 1-3 show the TRIAD solutions for the 
SAMPEX Euler attitudes. These attitude time histories are assumed to be the “true” 
attitudes for this test case. 

^e MME model error term d(t) is added only to the last three (angular momentum) 
states of the dynamic model, represented by Equation (9). The first four states (i.e. 
t e quaternions) are assumed to be perfect kinematic relationships so that no model 
correction is added to these states. This formulation avoids any difficulties encountered by 
the normalization constraint of the quaternions [2]. Therefore, no pre-conditioning of the 
estimator model for the normalization constraint of the quaternion states is required This 

formulation has clear advantages over the Kalman filter method for attitude estimation 
(see Reference [3]). 

Results indicate that the attitude time estimates given by the MME estimator are 
exact y identical to the TRIAD solutions, shown in Figures 1-3. Therefore, the MME 

aUUUde 6StimateS in thiS n ° n ' event case the estimates 
parallel the TRIAD solution throughout the entire time interval). 

The angular momentum estimates from the MME are shown in Figures 4-6 These 
trajectones are used to determine the instantaneous spacecraft angular rates resolved along 

^° d / pr ° pagate the fi uaternion s. The associated model error trajectories 

(<*(*)) from the MME estimator are shown in Figure 7. It is important to note that the 
correction is only applied to the angular momentum states. This formulation provides 
accurate MME state estimates of the Euler angles (see Figures 1-3), and also maintains 
the quaternion normalization constraint. 

The model error trajectories can now be used to correlate a linear or nonlinear 
correction to the SAMPEX dynamic model. To identify mathematical expressions that 
describe these trajectones, the Least-Squares-Correlation (LSC) algorithm is used [191. 

is algonthm develops a set of mathematical expressions that describe the model error 
histones as a combination of the state estimates. This algorithm can be implemented in 
wo ways: (1) the code can be allowed to form combinations of simple mathematical 
functions, or (2) a library of functions may be supplied by the user to augment the 
searc process (i.e., by supplying known functions from intuitive implementations or 
past studies). This library is a list of functions which the user expects will appear in the 
system under investigation. These functions may be common occurring functions from 
initial runs of the LSC algorithm. 

A method for identifying possible library functions involves plotting the model error 
trajectones versus the state estimates. These plots may offer significant mathematical 
insight on how to formulate library expressions. Plots of the second model error versus the 
rst and second state estimates are shown in Figures 8-9. From these figures, a possible 
functional form may be a Lemniscate geometric function with internal oscillations This 
geometric function is implemented into the library set (along with previous internal 
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functions). Table 1 contains example candidate functions obtained to this point The 
second model error candidate expression shows a high correlation coefficient o 0 . . 

Hgum 10 shows the second model error and the example candidate function that describe 
it g From the high correlation value, and the presence of only minor discrepancies in 
11” second model error is assumed to be correctly identified by the function shown 
in S Table 1. Note that this model error candidate expression has been identified without 

the use of any library terms. 

The next step in the study is to create a library starling with known attitude | dyn^l 
model components and external disturbances (e.g., aerodynamic torque, radiation torque 
orbit maneuvering torques, solar radiation pressure, etc). These can be used to further an 
more accurately identify nonlinear teims for the remaining model etror 
to obtain higher correlations). Once dynamical error models are obtained, they can be 
used to determine actual trajectories during anomalous penods such as Sun occultatio 
and/or measurement vector co-alignment. 



Figure 1 TRIAD and MME Euler 1 Angle Solutions 
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Figure 7 MME Angular Momentum Model Error Trajectories 


507 








Error 


Candidate Function 


Correlation 


-2.416 x lQ- 6 e fol a efo] 4 -2,120 x lQ-V^V^P +2.130 x 10~ 3 0.8647 

-4.581 x lQ-W’gj - 9.737 x 10~ 3 g f + 5.780 x 10~ 6 0.9963 

5.250 x 10~ 4 gjlf -5.163 x lO -5 # 3 ^^ + 5.163 x 10~ 5 0.8074 










.Ml 111 llill II 


SUMMARY 

A technique has been described which leads to algorithms capable of accurate attitude 
estimation in the presence of significant model error and/or 

satellites such as SAMPEX, attitude rate measurements are not available either by desig 
or by failure of existing rate gyros. The absence of rate measurements increases the 
estimation sensitivity to modeling uncertainty and measurement noise m the remaining 
relatively sparse attitude measurements. The technique described directly addresses the 
problem of attitude estimation without rate gyro data. 

Results using actual SAMPEX data and corrected models indicate that the technique 
described in this paper produces accurate estimates for both the spacecraft’s position 
and attitude rate. Also, the formulation described in this paper avoids any difficulties 
encountered when using quaternions to represent the attitude of *e satellite The 
new technique may be used directly as an estimator, or, as described in the paper, 
as a robust method of automatically obtaining accurate dynamic models for exisung 
satellites. Nonlinear candidate functions have been identified with fairly high correlation 
for SAMPEX. Later studies will utilize more candidate functions in order to obtain unity 
correlations for all model error trajectories. Therefore, these identified functions can 
used to propagate the model accurately in order to determine attitude and rate mouon 

during anomalous conditions. 
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